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NONLOCAL  CONTINUUM  THEORY 
FOR  DISLOCATION  AND  FRACTURE 

AoC.  Eringen 


ABSTRACT  -  By  means  of  linear  theory  of  nonlocal  by  means  of  classical  elasticity  (linear  or  ncn- 
elasticity,  solutions  are  given  for  some  Voltera  linear,.  Among  these,  we  mention  the  state  cf 

dislocations  situated  in  an  elastic  solid.  The  stress  at  a  sharp  crack  tip,  surface  ter.sicr. , 

stress  fields  are  determined  for  screw  and  edge  atomic  inclusions,  defects,  granular,  porcus  ar.c 

dislocations.  The  stresses  and  elastic  energy  composite  solids. 

are  devoid  of  usual  singularities  predicted  by  In  several  previous  papers  [cf.,  l-4],wchs-c 

the  classical  (local)  elasticity.  A  theory  is  shown  that  the  stress  fields,  due  tc  dislocations 

developed  for  continuous  distributions  of  dislo-  and  cracks,  predicted  by  the  nonlocal  elasticity 

cations  on  the  basis  of  nonlocal  incompatible  contains  n£  singularities.  In  fact,  tney  vanish 

elasticity.  Stress  fields  are  given  for  volume,  at  the  center  of  dislocations  and  at  the  crack 

surface  and  line  distributions  of  dislocations.  tip.  Moreover,  the  maximum  stress  occurs  at  a 

Celebrated  Peach-Koehler  formula  is  modified  to  short  distance  away  from  these  points.  By  equal - 

include  nonlocal  Green's  functions.  The  stress  ing  the  maximum,  stress  to  the  cohesive  stress 

fields  for  three-  and  two-dimensional  cases  and  that  holds  atomic  bonds  together,  a  physically 

for  the  anti-plane  strain  are  determined  for  line  realistic  fracture  criteria  was  estah lished- > 1 . 

distributions.  Calculations  are  carried  out  for  Ir.  the  classical  limit,  the  celebrated  Griff itr. 

the  uniform  distributions  of  edge  and  screw  dis-  criterion  is  obtained  with  the  dividend  that  the 
locations  along  a  straight  line  segment.  By  means  Griffith  constant  is  determined  without  any  acci- 

of  the  maximum  stress  hypothesis,  a  fracture  cn-  tional  assumption  on  the  surface  energy  which 

tena  is  introduced.  Calculated  theoretical  could  not  be  measured  to  within  any  reasonable 

strengths  are  in  good  agreement  with  those  based  accuracy.  Estimated  errors  ir.  such  measure-er.t* 

on  the  atomic  models.  Reduction  of  material  is  known  to  be  not  less  than  several  hundred  per- 

strength  with  the  presence  of  dislocation  line  cent. 

and  the  maximum  number  of  dislocations  are  given.  The  dispersion  curves  predicted  by  the  ncnlc- 

/V  cal  theory  are  nearly  in  coincident  with  those 

IT  IS  A  KELL- KNOWN  FACT  that  stress  fields  <  \  based  on  the  atomic  lattice  dynamics  and  obserx a- 

due  to  Volterra  dislocations,  contain  singulari-  tions  on  phonon  dispersions  [cf.  J-c], 

ties  at  the  center  of  the  dislocation  50  that  in  There  exist  ample  evidence  that  recer.tlv  de- 

a  small  region  around  the  center  (the  core  region1  veloped  theory  of  nonlocal  elasticity*'"*1  is  a 
classical  elasticity  fails  to  apply.  The  radius  proper  mathematical  model  which  can  eliminate  tr.e 
of  this  region  is  estimated,  usually  by  means  of  foregoing  difficulties,  extending  the  domain  of 

atomic  models.  Because  of  these  stress  and  energy  applicability  of  continuum  mechanics  to  physical 

singularities,  several  other  methods  have  been  phenomenon  with  internal  characteristic  lengths 

devised  for  the  prediction  of  fracture.  Clearly  at  the  molecular,  atomic  or  microstructura.sca.es. 

such  singularities  are  non-physical  and  a  proper  Based  on  these  observations,  we  expect  tnat  a 

model  should  eliminate  them.  nonlocal  theory  may  bear  fruit  m  dealing  se 1 

Phonon  dispersion  experiments  have  shown  that  stress  and  energies  of  dislocation,  loop?  fen  “~t- 

the  phase  velocity  of  plane  waxes  in  single  cry-  the  discrete  and  continuous  dist  nrutior.s  of  d:«- 

stals  depends  on  the  wa\e  length  so  that  disper-  locations.  Moreoxer,  the  onset  of  fracturt  a n_ 

sion  is  the  rule  rather  than  tiie  exception.  Vet  the  theoretical  strength  cf  solids  -a-  be  esti-.,- 

classical  elasticity  predicts  constant  phase  veio-  ted  by  means  of  the  continuum  theory  which  pe—;t 
cities  for  plane  waxes,  independent  of  frequency  extensions  tc  anorpno-s  solids  an.c  composites.  7:t 

and  wave  length.  rexsev  c  T:r,  of  the  prese-t  paper  ste-s  fre- 

There  are  many  other  physical  phenomena  in  tnese  observations, 

the  microscopic  scale,  that  cannot  be  precicted 


In  section  3,  I  present  a  summary  of  the 
linear  theory  of  nonlocal  isotropic  elastic  solids. 
In  section  3,  the  stress  field  and  energy  due  to  a 
screw  dislocation  are  calculated.  Both  turned  out 
to  be  devoid  of  singularities.  The  stress  field 
due  to  an  edge  dislocation  is  treated  in  section  4, 
In  section  S,  I  develop  a  theory  for  the  continuous 
distribution  of  dislocations.  Field  equations  are 
obtained  for  stress  functions  for  two  and  three- 
dimensional  state  of  strain. 

Green's  functions  for  the  infinite  solids  are 
obtained  in  section  b  leading  to  a  generaliced 
Peach-Koehler  formula  for  the  stress  field  due  to 
line  distribution  of  dislocations.  Results  are 
gratifying  in  that  they  contain  no  singularities 
so  that  calculations  can  be  carried  out  with  uni¬ 
formity,  for  surface  and  line  distributions  of 
dislocations  throughout  core  regions.  In  fact, 
stress  fields  due  to  a  uniform  distribution  of 
screws  along  a  line  segment  verifies  our  predic¬ 
tions.  The  theoretical  strength  of  a  single 
crystal  predicted  by  the  nonlocal  theory  is  in 
agreement  with  that  known  in  atomic  theory.  Also 
given  are  the  shear  strengths  reduction  due  to 
line  distributions  and  the  maximum  allowable  num¬ 
ber  of  screws  within  a  straight  line  segment  of 
length  2>  . 

2.  BASIC  EQUATIONS 


ar.  influence  function  pip.  c . 

which  brings  influences  cf  strains  at  •.aric.-s 
points  x'  to  x  ,  ir.  different  proportions . 
According  to  the'uttenuu-t -.c tc-.c hypo¬ 
thesis^,  it  assumes  the  maximum  value  at 
x’*x  ,  sharply  decreasing  with  the  distance  from 
x.~  From  Eq.  (3.3,,  it  is  clear  that  a  depenu: 
on  a  length  scale  e  .  Tms  is  an  internal  charac¬ 
teristic  length  which  may  be  selected  to  be  pro¬ 
portional  to  the  lattice  parameter  a  fcr  single 
crystals,  i.e. 

(3.5)  e  =  ec  a  , 

average  granular  distance  for  amorphous  bodies 
and  the  average  distance  for  fiber  composites, 
etc.  In  Eq.  (3.5.,  is  a  non-dimer.sitr.al  con¬ 

stant  which  can  be  determined  by  one  experiment. 

When  e  -  0  ,  Eq.  ;,5.3  must  revert  to  . 

This  implies  that  a  is  a  Ulnae  delta  secuer.cc. 
Thus,  in  this  limit  nonlocal  theory  reverts  tc  tr.e 
classical  elasticity  theory.  By  discreticir.g  Lc. 
(3.3),  it  can  be  shown  that  equations  of  nor.Iccal 
elasticity  also  reduce  to  those  of  atomic  latti.e 
dynamicsll.  . 

In  several  previous  papers’"’"’1'3,  I  gave 
special  representations  fcr  a  x'-.\  which  lea. 
to  excellent  predictions  ir,  accord  with  the  at:-, 
lattice  dynamics.  For  example,  fcr  the  twc-ci-e: - 
sional  case,  an  appropriate  kernel  is 


From  the  atomic  theory  of  lattice  dynamics 
and  experimental  observations  on  phonon  disper¬ 
sions,  it  is  veli-krowr.  that  the  stress  at  a 
material  point  x  in  a  body  depends  not  only  on 
strains  at  x  but  also  on  strains  at  al I  other 
points  x'  of  the  body.  In  linear  theory  of 
nonlocal'elasticity,  this  is  expressed  by  an  in¬ 
tegral  constitutive  equation  of  the  form  [cf.  3,4], 


(3.1.,  tv 


c,  .  (  x'  -x  :  e _ ( x * )  dv  tx  ’ ; 

k'-mn  -  -  mn  - 


where  c^  ^  are  the  material  property  functions 

which  depend  on  the  vector  x'-x  and  is 

the  linear  strain  measure  defined  by 


(3.3;  eu(x’;  •  4  [  ;x 


1 


wnere  u^yx.t,  is  the  displacement  vector.  For 
homogeneous  an.  isotropic  solids,  y 3 . 1 )  may  be 
simplined  to 


(3.5.,  tk>CA. 


,1  x'-x  i  (x';  dv(x'.. 


c 


wnere  ",  is  the  classical  local  stress 
x  • 

tensor  gi’.ti.  O'  Hoot.e's  ia- : 


k. 


e.  x 


(3.6) 


a,,  x 


where  is  the  modified  Sessel's  function, 

two-dimensional  lattices,  Eq.  3 .  c .  provides 
excellent  match  between  acoustical  disperse:" 
curves,  based  on  the  nonlocal  elasticity, 
based  on  Bom-karmar  theory  of  atone  lat 
namics.  In  the  entire  Brilioum  cone,  th 
is  less  than  6).  It  is  also  interesting 
that  for  the  infinite  media-  a  is  the  G 
function  of  a  linear  differential  operate 
the  case  of  Eq.  (3.6.  ,  this  means  that  a 


c  tr.esc 

error 
n.cte 
■er. ' s 
Ir. 

1 1 1 sties 


(3.-) 


(I  - 


x'-x 


vanishing  at  infinity. 

Under  the  raid  assumptions  of  vanishing  ncr 
local  effects  for  the  body  forces  and  couples, 
the  momentum,  balance  laws  of  nonlocal  elastccit 
reduces  to  Cauchy's  laws 

=  0 


is 


ete  summatio 
over  the  range  of  indices,  he  also  abreiiate 
partial  differentiation  with  respect  tc  \ 
ar.  index  following  a  comma,  ar.d  use  a  super-ms 
dot  to  indicate  the  time  derivative,  e.g. 


(3.s; 

tk;,k  *  '  “■ 

(3.9; 

V  s  :.k 

where 

C  is  the  mass  density 

body  force  density. 

As 

usual,  repeated  indices 

are  tnc  usual  Lame  con- 
,  tiu  i.emel  x'-x 


k  ,k 


i  r.  u  1 1 1  w . . 
starts. 


anu 


15 


r.\. 


(-.13; 


The  field  equations  of  nonlocal  elasticity 
are  obtained  by  combining  (3.1,,  (2.3'  and  (2.8,'. 

(2.10;  -  ,  c,  -  (  x’-.x  i  e  ix 1  j  dad 

kunn  -  -  ran  -  k 

•  V 


c.  (  x’  -x  ';  u  ,  (x1 1  dv’ 
kxmn  -  -  m,nk  - 


♦  C(f.  -  u. )  =  0 


where  a  superposed  prime  denotes  dependence  on 

x1  ,  e.g.  u'  =  u  (xM  ,  dv1  »dvix').  In  de- 
-  mm-  - 

riving  (2.10),  we  used  the  identity 


3  C 

r  ,  ,  .  k..  liu i  ■ 

w —  (c,  .  e  l  *  -  - — ; -  e 

ktmn  mn  ox^  mn 


ex'  kxmn  ran 

ie ' 
mn 

♦  C.  .  i — : — 

k-rar.  rxj( 


and  the  Green-Gauss  theorem  to  convert  the  first 
term  to  the  surface  integral  over  3C  in  (2.10). 
In  (2.10;,  the  first  integral  represents  the  sur¬ 
face  stresses  (e.g.  surface  tension;.  Consequent¬ 
ly,  nonlocal  theory  accounts  for  the  aurfsce 
phyeies.  This  important  asset  of  nonlocal  theory 
is  not  included  in  classical  field  theories. 

The  integro-partial  differential  equations 
(2.101  replaces  Xavier's  equations  of  classical 
elasticity.  The  displacement  field  u  is  to  be 
determined  by  solving  (2.10;  under  appropriate 
initial  and  boundary  conditions.  Boundary  and 
initial  conditions  on  displacement  and  velocity 
fields  are  identical  to  those  of  classical  elas¬ 
ticity.  Boundary  conditions  or.  tractions  is 
based  on  the  true  szrsss  censor  t,  .  ,  not  on 


,  i.e. 


(2.11, 


t,  n. 
k-.  i. 


in,.  ■ 


'  (n  :  . 


are  prescribed  boundary  tractions. 


For  mixed  boundary-value  problems,  to  avoid  a 
possible  overspecification,  care  is  necessary 
along  the  boundary  of  the  two  surfaces  on  one  of 


which  u.  and  or.  the  other  t 


(n 


are  pre¬ 


scribed14. 

If  we  assume  that  the  nonlocal  kernel  a  satis¬ 
fies  (2.-;,  then  for  the  homogeneous  and  isotropic 
solids  of  infinite  extent,  integro-differential 
equations  (2. 1C.  can  be  replaced  by  sir.p-ularZ^ 
piTZ'uVced  partial  differential  equations.  This 
is  achieved  by  noting 


'  “k,U  .,kk 

♦  (l-s:7-;.  (cf^  -  ;u.  =  0 

In  the  static  case  and  vanishing  body  forces 
reduce  to  Xavier's  equations 

(2.14)  P-Du^.  ♦  -  u,jkk  =  C 

However,  note  that  the  stress  field  is  deter— 
by  solving  (2.12,. 

3.  SCREW  DISLOCATION 

A  screw  dislocation  is  introduced  by  cute 
a  solid  along  the  plane  x,  =  0  ,  x^  >.  0  and 

troducing  a  constant  displacement  discontinue 
b  (called  Burger's  vector  along  the  xj-dire 
of  the  rectangular  coordinates.  In  this  case 
displacement  field  has  only  single  component 
u.(Xj,x,,t;  which  is  determined  by  selling 

(3.1;  7 


=  0 


The  stress  field  is  obtained  by  solving  2.12 
where  C-.  and  c.„  are  the  only  nor.-var.isu 
components  of  the  lical  stress  tensor.  These 
u.  are  given  byls: 

b  , 

(3.2;  u.  =  —  v  ' 

t3-3’  :31  =  -  F7  51n-'  -32-f^''cs: 

where  (r,r,c.  are  the  cylindrical  coordmat 
i.e. 

(3.4)  x,  =  r  cos  i  ,  x.  =  r  sin  r,  x.~ 

To  obtain  the  stress  tensor  t,  .  ,  we  de 
mine  the  solution  of  (2.12,  under  tne  cor.diti 
that  t.  must  vanish  on  a  circular  cylmdri 
surface*  of  radius  r,  2  '  6  '  ^  2  as  r  ■ 

Such  a  solution  was  giver,  in  a  previous  paper 

(3.5,1  t. 


‘■te 


This  is  regular  for  all  0<_r<_“  .  It  is  val 
in  the  core  region.  In  fact,  for  r  -  C ,  t.. 

vanishes  so  that  the  hoop  stress  is  cere  at" 
the  center  of  the  dislocation,  a  result  which 
in  accord  with  the  physics  of  the  problem. 

Classical  elasticity  solution  results  in 
case  of  E  = 0  ,  leading  to  a  stress  singular 
r*1  .  The  strain  energ'  per  unit  lengf.  of 
in  a  region  bounded  by  concentric  cylinders  o 
radius  r.  ,  and  R  ,  is  giver.  b> 
u  R 


.12 


;  1  -  c  * 


1-1.0 


:  i 


*  (xn-.R  r. 


x .  r.  C  -  x 


Contrary  to  the  classical  result,  this  has  no 
singularity  for  rQ  1  0  .  In  fact. 


is.-! 


IV0 


=  ^  [in ;R/2ej  ♦  K0 CR/e; ] 


which  makes  sense  on  physical  grounds.  For  large 
R,  Eq.  (3.7)  may  be  approximated  by 

(5.S,  I/L|  =  [in (R/2ej 


!V° 


*  (TC/2R)'4  exp(-R/£)] 


which  exhibits  the  dependence  of  the  energy  on 
the  sice  of  the  solid. 


(•J.s; 


I  f  k 


K  (f.  1 


C  ■  r/c 


Here,  I  (r)  and  k  (c)  are  modified  Bessel 
n  n 

functions.  By  superposition,  we  obtain 
(4.9)  0  =  i  Be'1  fjU,  ,elr  -  e'1'  , 


(4.10i 


i  Bc'^f,  if  e'"'»  f-(: .  e' 


4.  EDGE  DISLOCATION 

A  straight  edge  dislocation  in  a  solid  causes 
plane  strain  with  displacement  vectors  character- 
iced  by 

(4.1)  Uj  *  UjUj.x,);  u,  »  u,(Xj,.\,; ,  u.»0 

where  u^  are  the  rectangular  components  of  the 
displacement  vector.  The  displacement  component 
U’  undergoes  a  constant  lump  discontinuity  b 
along  the  half  plane  0,  x.,  =  0  .  The  classi¬ 

cal  elasticity  solution  of  this'problem  is  well- 
known^.  por  convenience,  we  write  the  classical 
stress  in  the  form 


(4.2) 

where 

(4.3i 


0.  =  Cjj  -  c,,  *  iBr'1 (e1  -  e'1')  , 

4_  *  °22'Cll*‘i  °i2  “  iBr"1(e"i'*e*';>i') 


B  =  ub/27(l-v) 


The  integration  can  be  carried  out  fcr  r.  =  l 
that 


(4.11) 


fj(C  =  c~l  -  k. 


From  Eq.  (2.12),  it  follows  that  the  true  stress 
field  satisfies 


(4.12)  f.(C)  =  |  I.Cf.  k . .  df 


*  !  K.(C-)  Klf.df 


The  stress  field  follows  fro-  .4.5 
uh 


ji,  -  sir.  -  *  r .  :  sir.  c? ,  , 


ii  4-(i -w  ,£  1 


(4.13)  t,,  .  -  r^-  [f,  ;  sir.  f-  f.  sir.  3- 


3,- 


ube 


1:  4- 


“y: .  —  ifi ;  co* r  *  fc 


(4.4)  (1  - e'V"){0,4}  =  {0,  ,tj 

where  6  and  $  have  the  forms 


*11**22  ’ 


*  “  t22'tll*2itl2 


In  polar  coordinates,  components  of  the  stress 
tensor  follows  from 


(J‘1J)  *rr**7f‘-2*  ^•trr*;i:, 


Thus,  we  must  find  the  general  solution  of 
2_2  -1  in1 

t4.b,  (1-c  .  jF  =  (£r  e  ,  n*:l,-j 


The  solution  of  (4,6)  which  is  regular  at  r=0 
and  r*x  is  found  tc  be 


Hence , 


(4.15; 


t  =  .  - — _  [f  - 

rr  4r  lV“  ‘f  -  -* 

t::  -  -  T.-T.—  .  [3f- .c  -  ]fir 


F  =  fnic;  e 


where 


t  f  *  *  f. 

rr  4-.  ]-.  c  L‘  1  *5 


Again,  we  notice  that  the  stress  f:clc  \  ar.isr.e? 
r  =  C  so  that  contrar'  tc  t:u-  classic-I  result. 


Erir.cer 


r 


I 


! 


! 


stress  singularity  is  present  at  the  center  of 
the  dislocation. 

It  should  also  be  noted  that  (a. IS  do  not 
reduce  to  the  formulas  of  the  classical  theory 
by  setting  £'0.  Tins  situation  is,  of  course, 
well-known  for  singularly  perturbed  differential 
equations. 


5.  CONTINUOUS  DISTRIBUTION  OF  DISLOCATIONS 

A  small  neighborhood  n(x)  of  x  in  a 
distorted  body  V  ,  may  be  relaxed  to  a  small 
neighborhood  N(X,i  of  the  image  X  of  x  ,  in 
an  undistorted  (.or  natural )  configuration'  V  , 
by  releasing  constraints  exerted  to  n(x)  by  the 
rest  of  the  body.  A  line  element  dx  at  x£n(x) 
can  be  expressed  in  terms  of  its  image  dX's  X  (xj 
by 

(5.1)  dx  *  A  dX 


The  strain  incompatibility  is  expressed  by 


(3.10, 

where 
and  it 

(3.11) 


lik  .,mr.  in  ,;m  k.. 

is  called  the  ino:~~.  cr i\i~. 
is  given  by 

n,  .  »  i  (c,  a  .  ♦  c .  a  . 

kv  -  krar.  n-.,m  vrar.  nx.r. 


All  of  these  results  are  well-known  in  classical 
theory  (cf.  [lb]  . 

Nonlocal  theory  stipulates  that  the  stress 
strain  relations  is  giver,  by  (Z.  1..  For  homoge¬ 
neous,  isotropic  solids,  under  a  mild  assumption 
on  the  nature  of  the  attenuating  kernel,  we  ha.  e 
(2.12;  ,  i.e. 

(3.  Id)  (1  -E'Vjt  .  =  /.  e  .  -  2d  e.  . 

k  a  rr  n '  f ■ 


where  A(X,  is  called  the  elastic  distortion. 
he  assume'that  A(X.  is  continuously  differen¬ 
tiable  and  possesses  a  unique  inverse,  so  that 
-1 

(5.2)  dX  =  A  dx 


From  this  we  solve  for  e^ 


where  .  =  '•  is  Poisson's  ratio.  If 

substitute  Lq.  (5.15  into  .3.1..,  we  ottair. 


Consider  a  smooth  surface  S  in  U  bounded 
by  a  closed  curve  C  .  The  true  Burger's  vector 
b  of  the  dislocations  piercing  through  S  is 
3efined  bv 

f  -1 

(5.3)  b  *  0  dX  =  0  A  dx  *  j  a  n  da 
C  ~  C  S 

where  n  is  the  unit  normal  to  S  ,  the  posi¬ 
tive  sense  of  C  being  counter-clockwise,  when 
sighting  along  n  .  Mere,  a  is  called  the 
true  dislocation' density 


(5. Id.  U-c*7‘ 


t,  ■  *  —  it. 
k  ■  1  -  rr,> 


The se  equations  m..st  be  solved  for  t-  ur.ic: 
are  sufciect  tc  tne  equilibriur  conditio-. 

V,k  - 

Following  kroner's  classical  aprr;*;:.*  ,  re j:  r 


ing  tne  Beltrami  solution  : f  3.25  , 

«  t1  t  Sf  T 

-1 

-l 

(5.4  ' 

a 

*  curl  A  , 

a_ik  "  €kmn  A.in,m 

i5.lt  t.  *  7“ . . 

K  K  - 

For  small 

distortions , 

we  can  write 

- 1 

1-.  ' ' rr.k  ’t 

T  k  • 

(3.5; 

\ 

1  *  :k:  *  ■‘k. 

•  '  :k:  '  ak. 

so  that 

sut-ect  tc 

(5.6, 

a;k  "kmr.  uim,n 

;  5 . 1  “  ^ .  *  C 

From  this. 

it  follows 

that 

Substituting  5.1b  into  3.1-  , 

w  e  ;  r  t  a ;  * 

(5.  *) 

\ik,k 

*  0 

(5. 15  ,1  -  e*T* 

K  >  t 

Tne 

linear  strain 

tensor  e.  and  rotation 

. 

tensor 

wk 

are  giver. 

by 

through  ,3.;:  ,  we  calcuiatt  the  mo 

c -p atic.  .it 

i.5.  S. 

i 

ek.  1  : 

(a  ♦  j  . 

^  ■  *.x 

tensor  .  Tne  scl-t:or.  cf  5 . 1 s 

K  - 

give  xk.  • 

vs...  tr.cr 

As  expectec.  L - .  ,5.1s  rec..;ec  t  c  tr-.  dj 
cal  e c  ja t 1  or. ,  * ;< e r. 


pert-rbec  partia.  cirrerer 
tain  t r.e  sc Ijti or.  ::  5...' 


::r.c  jr 


w.  . 

k  • 


ec .. 


function,  G,  ix.:  ,  which  must  satisfv 
K*.nn  -  ~ 

1,5.19.  ( 1  -  t  ;  ■  G.  *  1 1, x-£  '  - 1  ■  - „ 

k  a-  mn  -  -  k  *  mn 


<:  -  £"",g 


When  G,  is  known,  then  the  solution  of  (5.18) 

is  givdrs? 


(S.:o)  XU  *  }  ^(5)  d'(Ci 


In  the  case  of  the  plane  strain,  introducing 
Airy's  stress  function  by 

(3.21)  tu  »  *  ,  t,,  «  4>  n  ,  tj,  =  -  i  j. 


we  obtain  an  equation  replacing  (3.18! 

,,  ,2-:,  .  ,  „ 


(5.23) 

(1  -e:T; 

where 

(3.23,:  h 

=  n33  -  3 

(3.24)  a,. 

a  J 

=  a2l.2-! 

depend  or. 

Xj  and  x. 

In  the  case  of  the  anti-plane  strain,  equa¬ 
tions  of  equilibrium  are  satisfied  if 

(5>:5  li3  *  IxT  t:3'=-% 

and  we  obtain 


(5.36,  (1-C 

where 


*  u  a.- 
00 


a..  =  a.  ,  -  a.,  , 
oo  ol,-  o. ,  1 


For  the  infinite  space,  h  is  gi.e.n  by 
(6.3;  H  •  -  .\  *  ;  5~ 

In  spherical  coordinates,  using  the  opera 

d  ,  :  c 


dr  ‘  dr' 


we  obtain 


(6.4)  G('x-q|  j  B  4--”xI - —  exp  -  x-3  it, 


(b,5j  G(  x-c  ,  =  -  \  - "  5“  ,  £sC 

where  we  also  determined  instant  tc  re 
tkk  regular  at  j  =  i  , 

The  solution  of  (5.  for  the  ir.fi- 
media,  is  given  by 


(6.6-,  >k. 


which  satisfies  conditions  \5.1~  on  acco 
(5.*)  and  (5.11  ■.  Upon  substituting  fret 
this  mav  be  expressed  as 


(6."i  Xkx(xi  m  “  ciik  1  a.: ''  ~  dv 

'  V  ' 

♦  it..  ,  a  . ,  ( ’ "  — —  d \ 
2  !-■  1  ik  ~  rx 


6.  GREEN'S  FUNCTIONS:  STATE  OF  STRESS 


The  determination  of  the  stress  fields  aris¬ 
ing  from  the  continuous  distributions  of  disloca¬ 
tions  requires  that  we  obtain  the  Green  function 
for  the  stress  functions  .  ,  i  and  I  .  Here 
we  determine  Green's  functions  for  solids  of 
infinite  extent. 

(i)  THREE-DIMENSIONAL  INFINITE  SOLID  -  The 
operator  T*  is  invariant  under  the  rotations 
of  coordinates.  For  the  infinite  space,  we  seek 
a  solution  of  :5.1S  which  depends  upon  1 x-i 
on  1  y ,  i .  t . 


where  we  used  the  Green -Gauss  theorem  and 
surface  term,  at  infinity  to  :ero. 

Eq.  (6.")  may  be  used  to  study  varicu 
cases  involving  surface  and  line  distribu 
dislocations.  For  example,  for  a  line  di 
tion  of  dislocation  along  a  closed  curve 
obtain 


(b.S;  4-  £...  b  ®  As¬ 
h'  .ilk;  x 


(6.  Ij  (1  - c~7~ J  7  G  *  : ix-3 


Since  the  operators  (1  -e 
we  set 


and  7  commutes,  where  b.  is  the  Burger's  vector  per  un: 

of  C  aiid  ds  is  tne  element  of  arc. 

i 


The  stress  field  due  to  a  line  distribution 
of  dislocation  is  obtained  b>  substituting  'o .S 
into 

^.9  0  ^V 

c 

♦  if  ^i  -  ^.iV^r1 


This  result  is  identical  to  the  Peach-Koehler  3’ 
formula  with  the  modification  that  here  G  is  the 
nonlocal  Green's  function  i.b.4 )  with  c  *  u.  The 
most  interesting  new  feature  of  kb.9;  is  that  it 
does  not  exhibit  unbounded  stress  fields  and  ener¬ 
gies  at  any  point  on  C  or  elsewhere. 

(ii,  TWO-DIMENSIONAL  INFINITE  PLANE  -  Green's 
function  in  this  case  can  be  fqund  to  be  similar 
to  decomposition  (,b.L;  with  T*  and  H  ,  given  by 


r  dr  'T  dr" 


H  =  -  T7  In  x-C 


(b.i:,  G,  x-J  ;  *  yr  K0l.x-c;/c; 


(x-C/* 


The  stress  is  calculate!  b;-  usir, 


(b.lOj  t„,  =  -lu  _  (G  j . . b.dt .  -  G  ,  b. 


t  =  ;g  , . -b.d-: .  -  g  ...t 

1-  ;  f 1 i -  -  . 


where  indices  after 


>mr.a  denote  cartes,  cer;  . 


tives  with  respect  to  ,  e.g. 


(iii ,  ANTI-FLAN-  STRAIN  -  The  Greet 
for  the  differential  operator  :r.  i.it 
tained  to  be 


(b. 1 G i  x-i 


-  A  -  X  - :  '  C  ]  , 

l  -  -  ' 


X-, 


(b.lSi  G (  x - i  ,  =  -  Lr.  x -  ‘ 


The  stress  function  :  is  giver,  by 


lb. 1  Si  G(  x-i. ; 


(x-:  i*  (x-;/ 


in  ( :  x-5  i  > 


(6.19  :  iX 


G,  x-: 


The  s.ress  fielc  is  found  to  be 


where  K,(:,  is  the  modified  Bessel's  function. 

0 

Airv's  stress  fuention  is  obtained  to  be 


(6.:oi  t. 


t-i  a-:,  dl- 


io.  14  i(\ 


r3G  ... 
;  tir-  ar5U; 
i  i 


-  tr  ais!^Jda 


w  ,  -  jr  J- 


where  we  used  the  Green-Gauss  theorem  and  set  a 
line  integral  tc  :erc  at  infinity.  For  s  Zi'.c 
distribution  of  dislocations  in  the  x_*0-plane, 
we  will  hare 


For  a  line  distribution  of  dislocations  or 
x.  =  0-piane,  the  stress  field  is  giver,  by 


(-.15  i  x  »  -  b-'!  dr. 


rr-  o- u' 


lr.  plane  polar  coordinates  ir,~  ,  we  nave 


i  hr- 


b(i;  ds 


-r—  bC  i  ds 

-T  - 


7.  STRESS  DISTRIBUTIONS 

Here,  I  present  some  results  on  the  stress 
field  due  to  continuous  distributions  of  disloca¬ 
tions  along  a  line  segment. 

u:  EDGE  DISLOCATION  ALONG  A  LINE  SEGMENT  - 
The  stress  field  due  to  an  edge  dislocation  can 
be  calculated  by  using  (0.16'-.  In  this  case, 
b.*C  and  we  consider  b^  *  const,  distributed 

along  a  line  segment  x-  *  0  ,  Xj  <  ■  .  Using 
t.ne  fact  that  G  ,  -  3G  ri.  *  -  ?G, hi-  ,  the  inte¬ 
grations  m  (o.lty  are  performed  readily 


tK  =■  :-t]  [g  ;;ur  -  g  . 

t;:  =  :ubj  [g^jUj-  -  G(,,u;')  , 

ti:  *  *  ;'-bi  tG,i:'bi  *  c,i:(£:» 


\-  v  x" 


Tne  shear  stress  gi  er.  by- 
expressed  in  nor.-airer.sicr.al  : 


(-.5  T,  *  t,. 


('.6;  td  *  -b  I-  *  -b,.N 


(in  SINGLE  SCREts  -  For  a  single  screw, 
shear  stresses  can  oe  obtained  b>  surst :  t -:  :r  . 

v-'-'  b-;  *  bc  :  i 

into  ,c.;C  ,  where  is  the  Dirac-cc 1: u 

measure.  Tms  lead;  tc 


r  u  x.1  -  •  .  -  1  « 

t.  -  »  -fc  — - *■  [.  -  -  h-  r  .  , 

,  -. c  *  -  °  i  "  r* 

t..  *  -  .  =  tit”  — F-  [1  -  t  K.  r  :  ] 

-U  0  ;x,  .  -  C 


where  r  =  (xt  *  xt  :  .  Tms  result  is  ident: 

to  Eq.  ( 3. S ,  obtained  different!.',  in  polar 

dinates . 

For  a  single  screw,  the  shear  stress  Ec, 
may  be  expressed  in  non-dimensional  for" 


=  [  x,-L.-  -  a!]'*  e  , 


[(x.-L  *  *  xt ]'*.  £ 


Green's  function  G  is  giver,  by  (o.ld 


(ii  SC  REF.  DISLOCATION  -  For  a  uniform  distri¬ 

bution  of  screw  dislocations  along  a  line  segment 
x.  *  0  ,  Xj  <  1  ,  through  ib.21,,  we  find  that 

.-.3  t,.  =  -  b[G(:  ,  -  Gv;;i] 

where  G  is  given  by  (o.l",  and  : ..  and  by 


If  the  number  of  dislocations  is  N  over  a  dis¬ 
tance  i  ar.d  b(^  is  the  atomic  Burger's  vector, 
we  ha.e  for  the  macroscopic  Burger's  vector. 


0  *  b^N 


The  stress  field  given  by  Eq.  is  displa. 
graphically  in  Fig.  1.  It  has  no  singularity, 
fact,  contrary  to  the  prediction  made  ir.  c.ass 
elasticity,  TcU  vanishes  at  :  *  0  .  The  no, 
mum  stress  occurs  at  b  -  1.1  and  is  g:\er  by 


t  -  =  C,  3993  -r^- 

crma  x  _"E 


The  fracture  will  begin  when  t.^  «  t  =  tr.e 

cohesive  stress  that  holds  atomic  bonds  together. 
Therefore,  nonlocal  theor>  predicts  that 


Note  that  while  this  is  or.!'  ate-: 
for  single  crystals,  neiertheless , 


Erin  get. 


1c 


finite  distance  for  amorphous  ar.d  composite 
materials .  It  depends  on  the  internal  character¬ 
istic  length  a  (more  precisely,  e^a  . 

Ck  Frasture  Sri teric^.:  ."he  fra? r  errors 

exc-:  the  -airc-u-  shezr  stress  si  se\  c- 
\  1 '  rcccht s  the  jc Due  eg' the  roheetve 
stress  thdt  holds  atetir  bonds  together. 

Kith  this  criterion,  the  maximum  stress  hypothesis 
for  fracture  is  restored  for  microscopic  and  atomic 
phenomena  as  well. 

If  we  write  h*e/0.  3993.  Eq.  C~.ll)  agrees 
with  Frenkel's  estimate  of  tiie  theoretical  strength 
of  single  crystals,  based  on  atomic  considerations 
(cf.  kelly  [18],  p.  12'.  In  fact,  if  we  use 
£  *et,a  =  0.39a,  we  find  for  the  single  aluminum 
crystal , 


t  /•„  *  o.i: 

c 


{AC:  [Ill]  <  llo >  } 


This  is  very  close  to  the  theoretical  strength 
t  /•„  =  0.11  based  on  atomic  models. 

In  the  case  of  line  distribution  of  screws 
with  constant  Burgers  vector,  the  non-dimensional 
shear  stress  is  displayed  in  Fig.  2  for  various 
values  of  Y  .  Here  we  observe  that  the  shear 
stress  is  maximum  near  the  end  points  of  the  line 
segment  Xj  ”  1,  x,  *  0  .  It  is  located  slightly 
outside  of  the  end  pSints  and  for  y  >_3  it  is 
very  close  to  the  end  points  (see  Table  1.1.  Again, 
contrary  to  the  classical  elasticity  solution, 
there  is  no  singularity  at  x  j  *  x.  *  0  .  Be¬ 

havior  of  T;  is  governed  basically  by  the  first- 
term  in  C.Sj  except  near  x  =  ♦_  1  .  At  x  =  1  ,  we 
have 

.....  ub0N  in-.  »B 


If  we  write 


B  =  b0N 


h  =  cy  '  ;n  Y 


Eq.  (".13,  may  be  interpreted  in  terms  of  the  slip 
of  atomic  layers  of  distance  h  by  one  narroeaepir 
dislocation  of  Burger's  vector  B  .  The  ratio  of 
the  cohesive  stress  t^  for  the  line  distribution 
to  that  of  a  single  screw  is  given  by 

[Table  1:  Maximum.  Shear  Stress  and  its  Location  i 
for  a  Uniform  Distribution  of  Screw  | 
Dislocations  Along  a  Straight  Line 
Segment 


l.AAfc  1.19"  1.103  1.039  1.000  1.000 
0."-i“S  1.0501  1.5008  i.0651  D.30De  3.995“ 


This  gives  tl 
presence  of 
ly  along  a  si 

Since  t  ,  <  1 
d  — 

tions  is  give 


he  shear 
3.N  disl 
traight  1 


stress  re  a  - 
ocati cr.s  di 
me  segment 


cr.  cue  : : 
ibutec  ur 


Of  course,  this  n-jmber  will  have  to  be  mocifee. 
when  the  distribution  is  not  uniform. 

The  stress  fields  due  to  a  uniform  distnc. 
tion  of  screw  dislocations  along  a  circle  are 
given  in  another  publication. 
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